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conductivity in quantum non-degenerate collisional plasmas under arbitrary degree 
of degeneracy of ,he electron g a s L received. The Wigner-Vlasov-Boltzmal kinetic 
equation with collision integral in BGK (Bhatnagar, Gross and Krook) form in 
coordinate space is used. Various special cases are investigated. 
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INTRODUCTION 



In the present work formulas for calculation of electric conductivity and 
dielectric function in quantum non-degenerate collisional plasma under arbitrary 
temperature, i.e. under arbitrary degree of degeneration of the electron gas 
are deduced. 

During the derivation of the kinetic equation we generalize the approach, 
developed by Klimontovich and Silin |T]. 

Dielectric function in the collisionless quantum gaseous plasma was studied 
by many authors (see, for example, [[Q [TP]). 
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In the work [6] , where the one-dimensional case of the quantum plasma is 
investigated, the importance of derivation of dielectric function with use of 
the quantum kinetic equation with collision integral in the form of BGK - 
model (Bhatnagar, Gross, Krook) [TT], [12] was noted. 

The present work is devoted to the performance of this task. 

A dielectric function is one of the most significant characteristics of a 
plasma. This quantity is necessary for description of the skin effect [13], for 
analysis of surface plasmons [H], for description of the process of propagation 
and damping of the transverse plasma oscillations [10], the mechanism of 
electromagnetic waves penetration in plasma [9], and for analysis of other 
problems of plasma physics P2], [IB], \T7\, PS] and [19]. 

Kliewer and Fuchs were the first who have noticed pE], that the dielectric 
function for quantum plasma deduced by Lindhard in collisional case does not 
pass into dielectric function for classical plasma in the limit when Planck's 
constant h converges to zero. This means, that dielectric Lindhard's function 
does not take into account electron collisions correctly. Kliewer and Fuchs 
have corrected Lindhard's dielectric function so that it passed into classical 
one under condition h — > 0. 

In the works [H] , [TB] the dielectric function received by them was applied 
to consideration of various questions of metal optics. 

In the work [5] the correct account of collisions in framework of the 
relaxation model in electron momentum space for the case of longitudi- 
nal dielectric function has been carried out. At the same time the correct 
account of influence of collisions for transversal dielectric function has not 
been implemented till now. 

The aim of the present work is the elimination of this lacuna. 

1. KINETIC EQUATION FOR THE WIGNER FUNCTION 

Let's consider the Schrodinger equation written for a particle in an elect- 
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romagnetic field in terms of density matrix p 

ih^ = Hp-H*'p. (1.1) 

Here H is the Hamilton operator, H* is the complex conjugate operator to 
H 1 H*' is the complex conjugate operator to the H, which forces on primed 
spatial variables r'. 

Hamilton operator for the free particle which is in the field of the scalar 
potential U and in the field of vector potential A, has the following form: 

(P " -A) 2 

H = ^ + eU = 

2m 

= fi-2^< pA + Ap > + 2i A2 + eC7 - ™ 
Here p is the momentum operator, p = —ihV, e is the electron charge, 

m is the electron mass, c is the light velocity. 

Let's rewrite the Hamilton operator (1.2) in the explicit form: 

H = -—A + —(2AV + VA) + —,h 2 + eU. (1.3) 
2m 2mc V / 2mc z 

Complex conjugate to the H operator H* according to (1.3) has the form 
h 2 A ieh / . _ _ . \ e 2 



H* = A - 



- — 2AV + VA +- — ,A 2 + eU. 

2mc \ J 2mc z 



2m 2mc 
Hence we can write down for Hp: 

fz 2 icfi / \ 

Hp = -2^ Ap+ 2^ ( 2AV " + " VA ) + W A " + eUp ^ 
and for H*'p: 

H "? = -i A '? ~ S-c ( 2AV '" + " V ' A ) + 5i A ' 2 " + eU ' p - ^ 

Operators V and A from Eqs (1.4) and (1.5) force on unprimed spatial 
variables of the density matrix, i.e. V = Vr, A = Ar. In the operator H*' 
it is necessary to replace the operators V = Vr and A = Ar by operators 
V' = Vr' and A' = Ar, in addition we introduce the following designations 

A' = A(R/,f), U' = U(R',t). 
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Let's find the right-hand member of the equation (1.1), i.e. difference 
between relations (1.4) and (1.5): Hp — H*' p. According to (1.4) and (1.5) 
we have: 

h 



H »- H "o = -^-*») + 



2m 
ieh r 



2mc 

„2 



2^AVp + A'V'p) +p(vA + V'a) 



+ 



2mc 2 



A 2 (R, t) - A 2 (R/, t) + e[U(K, t) - U{R', t)]p. 

The connection between density matrix p(r, r', t) and Wigner function 
20] . [ 121] . |22] f(r, p, t) is given by the inversion and direct Fourier conversions 



/(r, M = / P(r + f , r - |, t)e-^ 3 a, 



p(R, R', t) = 



1 



;27r^) : 



/( '?l±^ 5 p 5 A^p(R-R')/^ 



The Wigner function is analogue of distribution function for quantum 
systems. It is widely used in the diversified physics questions. Substituting 
the representation of the density matrix in terms of the Wigner function (1.2) 
into the equation for the density matrix (1.1), we obtain 



lH di = H 



1 



[27rhy 



-H 



*/ 



5l±^,p',^ e V(R-R')/^3 p / 



Let's use the equalities written above. Thus the right-hand member of the 
previous equation we may present in explicit form. As a result we receive the 
following equation: 

1 f 



.,dp 
in 



dt {2nhf 
ieh r 



ih .„ „ ieh 
— PV/ + — - 

m zmc 



divA(R, t) + divA(R',t) 



f- 



+ 



2mc 



A(R, t)+A(R',t) 



mc 



A(R,*)-A(R',*) p'/ + 
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+ 



2mc 2 



[A\R,t) - A 2 (R', /;)]/ + e[U(R,t) - U(R',t)]fj e ip ' {n - R ' )/h d 3 p' . 

(1.6) 



In the equation (1.6) we will put 

R = r + -, 

Then in this equation we obtain 
R + R' , 



R' = r - -. 

2 



/(ii±^,p',t) e V(H-R')A = /(r , p ^) e Va/, 

Let's multiply the equation (1.6) by e~ ipa ~/ h and let's integrate it by a. 
Then we will divide both parts of the equation by ik. As a result we receive 



df 

dt 



m zmc 



A(r + |,*) + A(r-|,*) V/+ 



mch 



A(r + |,t)-A(r-|,t) p'/+ 



2mc 



cl cl 

div A(r + -, t) + div A(r - -, t) 



It 



2mc 2 h 



A 2 (r + |,t)-A 2 (r-f,t) 



/- 
/- 



"ft 



f/(r + |,t)-f/(r-|,0 



/ e 



^(p'-p)a/ft^_^V 
(27T/i) 3 ' 



(1.7) 



On the left-hand side of the equation (1.7) we have / = /(r, p, t), in the 
integral we have / = /(r, p', t). 
We consider the integral 



p'(V/)e 



73^/ 



p' ,, a i, d :i «<I"P = v / / p'/e^p'-P*//' 



(2tt/i) J 



d 3 a 

(27T^) 3 



= p'Mp'-p^p^pV/^p). 
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Two following equalities can be verified similarly: 

e 



A(r,t)(Vf(r,p\t))e'V-^i^ = 

mc {ZTTh) 6 

= — A(r,*)V/(r,p,*), 

mc 



and 



= — (divA(r,t))/(r,p,t). 

mc 

Then the equation (1.6) can be rewritten as following 

% + - (p - - A) V/ - — (div A(r, *))/(r, p, t) = W[f}. (1.8) 

In the equation (1.8) the symbol W[f] is the Wigner — Vlasov integral, 
defined by the equality 



W[f] = 




2mc 



A(r + |, *) + A(r - |, *) - 2A(r, *) V/ + 



ic 



mch 



A(r + |,*)-A(r-|,*) p'/+ 



2mc 



divA(r + ^,t) + divA(r- ^,£) -2divA(r,£) 



/- 



2mc 2 h 



A 2 (r + |,t)-A 2 (r-|,«) 



/- 



tf(r + |,t)-tf(r-|,t) 



«(p'-p)a//i 



d 3 a dV 

(27T/l) 3 ' 



The energy of the particle is equal to 



£ = 



J_( p _£ A ) 2 + e [/. 



Then the velocity of the particle v is equal to 

d£ If c . 
v = — = — p--A 
op m\ c 



(1.9) 
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besides, 



Vv 



div A. 



mc 



Hence, the left-hand part of the equation (1.9) equals to: 

% + - (p - -A) V/ - /— div A = % + vV/ + / VA = 

ot m\ c / mc ar 

-g + V(v/). 



Therefore the equation (1.9) can be rewritten in standard for transport 
theory form 

g + V(v/) = W\f}. (1.10) 

In the case of collisional plasma we may write the kinetic equation (1.10) 
as following 

^ + V(v/) = 5[/,/] + W[/]. (1.11) 
In the equation (1.11) the symbol B[f, f] represents the collision integral. 

2. RELAXATION MODEL OF KINRTIC EQUATION 



Under the electron scattering on impurity we will consider the equation 
(1.11) with collision integral in the form of relaxation r-model [TT], [12] : 



f + V(v/) 



/<"> - / 



r 



+ W[f\. 



;2.i) 



In the equation (2.1) r is the mean time between two consecutive collisions, 
r = 1/^, v is the collision frequency, is the local equilibrium Fermi - 
Dirac distribution function, 

*£ - "M-i 



/ 



(o) _ 



1 + exp 



M 



Here /c^ is the Boltzmann constant, T is the plasma temperature, £ is the 
electron energy, /i is the chemical potential of electron gas. 

In an explicit form the local equilibrium distribution function has the 
following form 

|2 1 -1 



/ (0) M) = 



1 + exp 



[p-(e/c)A(r,t)_L + eC/(r,t)-/i 



2mkBT 



k B T 
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We introduce the dimensionless electron velocity C(r, £), scalar potential 
4>{r 1 i) and chemical potential a: 

, v(r, t) , . . eU(r.t) ll 



1 777/ 

where i>t = —j= is the thermal electron velocity, (3 = 



Vf3 J ' ^ 2fc B T' 

Now local equilibrium function can be presented in terms of the electron 



velocity as follows 

/<°> (r, t) = [l + exp f + n - 

J v ' ' L F V 2/^T k B T J 



or, in dimensionless parameters, 

1 + exp [C 2 (r, t) + 0(r, t) - a\ 
We designate x = ol — 4>- Then we have = 



The quantity x is defined from the conservation law of number of particles 



fdn F = J f {0) dtt F . 



Here dQ F is the quantum measure for electrons, 

Let's note, that in the case of constant potentials U = const, A = const 
the equilibrium distribution function (2.2) is the solution of the equation 
(2.1). 

Let's find the electron concentration (numerical density) N and mean 
electron velocity u in an equilibrium state. These macroparameters are defined 
as follows: 

N(T,t) = J f(T,p,t)d£l Fl 

«M) = N * r ^ J v(r,i)f(r,p,i)dn F - 
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For calculation of these macroparameters in equilibrium condition it is 
necessary to put / = f(°\ where /(°) is defined by equality (1.13). We 
designate these macroparameters in equilibrium condition through N(°\r, t) 
and u(°)(r,t). 

Let's carry out the replacement of the integration variable p— (e/c)A(r, t) = 
p' in the previous equalities. Then, passing to integration in spherical coor- 
dinates, for numerical density in an equilibrium state we get: 

JV<°> = ^|/ 2 (a - 4>), (2.3) 



where 



2 

tic dtjc 



/x doc f 
T~. , 2 : r = / x 2 f F (a - (j)) dx. 
1 + exp(x z + (p - a) J 
o o 

In the same way, as for numerical density, for mean velocity in equilibrium 

state we derive 

u (0) M) = ^y J v(r,t)fW(r,p,t)<m F , 
or, in explicit form, 

u<°>(r,t) = 



(0)/ - I [p-(e/c)A] d 3 p 



N(°)(27Th) 3 J 1 r(p-(e/c)A) 2 eU - /i 



1 + exp 



+ 



After the same change of variables p — (e/c)A(r, t) = p' we receive: 



u (»)( rt) = I / til -o (24) 

U " JVC) (2^)3 J _ r / eC /_ M - ' 



1 + exp 



+ 



-2kBTm 

So, the electron velocity in an equilibrium state according to (2.4) is equal 
to zero. 

Let's note, that numerical electron density and their mean velocity satisfy 
the usual continuity equation: 

dN 



dt 



+ div(iVu) = 0. (2.5) 
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For the derivation of the continuity equation (2.5) it is necessary to integrate 
the kinetic equation (2.1) by quantum measure for electrons dftp and to use 
the definition of numerical density and mean velocity. Then it is necessary 
to use the conservation law of number of particles and to check up, if the 
integral by quantum measure dQ p of Wigner — Vlasov integral is equal to 
zero. Indeed, we have 

/ W[f] ^k = 2 // { " ' } eiPW ^ (a) Mp ' = 



= 2 




d 3 p - 0, 



a=0 



as after some algebra, 



{■} 



= 0. 



a=0 



Here the symbol {• • • } means the same expression, as in the right-hand 
side of the equation (1.9). 

Let's note, that the left-hand side of the kinetic equation (1.11) or (2.1) 
takes standard form for transport theory under the following gauge condition: 

divA(r,f) = 0. (2.6) 

Thus, i.e. in case of gauge (2.6), the kinetic equation has the following 
form: 

df 



at 



+ vVf = B[fJ] + W[f}. 



(2.7) 



Here the Wigner — Vlasov integral equals to: 



W[f] = 




2mc 



A(r + |, t) + A(r - |, t) - 2A(r, t) Vf + 



le 



mch 



A(r + |,t)-A(r-|, t )]p7-£ 5sL 



A 2 (r+|,t)-A 2 (r-|,t) /- 



Ti 



U{T+\,t)-U{T-lt) 



i(p'-p)a/fi 



d 3 a d 3 p' 
{2irhf ' 



(2.8) 
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3. LINEARIZATION OF THE KINETIC EQUATION AND ITS 

SOLUTION 



Let's consider the kinetic equation with collision integral in the form of 
r-model and suppose, that the scalar potential is equal to zero: U (r, t) = 0. 

We take vector potential which is orthogonal to the direction of the wave 
vector k : kA = in the form of a running harmonic wave: 

A(r,f) = A e^ kr -^. 

We suppose that the vector potential is small enough. This assumption 
allows us to linearize the equation and to neglect terms quadratic in electric 
field. 

Then the equation (2.7) can be reduced to: 



(3.1) 



In this case chemical potential is equal to a constant. 
In the equation (3.1) local equilibrium Fermi - - Dirac distribution is 
simplified as following: 

/ (0) = [l + exp(c 2 (r,f)-a X ' 



(3.2) 

The Wigner — Vlasov integral (2.8) also can be simplified essentially and 
has the following form: 



W[f] = 



le 



mch 
We notice, that 



a 



A(r + = t) - A(r - ^t)] P >fe^-^^0. (3.3) 



A(r + |,t)-A(r-|,i) = A(r,i) 



Jka/2 _ -ika/2 



Calculating the integral in (3.3), we find, that 



le 



ika/2 _ „-ika/ 2 



,«(p'-p)a, 
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The internal integral is equal to: 
1 





kft- 






kfti 


| exp yi 


p'-p + Y. 


ft/ 


— exp ^ 


P'-P + T . 



^(p'- P+ f)-<p'-p-f) 

We calculate the Wigner — Vlasov integral 



ftk\ 



= A(r,*) 



mcft 



W[/,A] = 

<5(p' - p + y) - 4 (p' - p - y) P'/< r > p' *) d V 



= A(r ' ^ i ( p " t) /(r ' P ~T'^~ ( p + t) /(r ' P + T'^ 



hk 



ftk^ 



hk 



A( t) ^ C { 



<■/ . „. ftk ' 

/(r,p-— ,*)-/(r,p + — ,*) 



ftk 

T 



ftk 



ftk "ii 

/(r,p-^,t) + /(r,p + — ,t)jj 



where 



= A <''0('+-'-)' 



f± = /(r,PTy,i). 



Consequently, the Wigner — Vlasov integral is equal to 



le 



mcft 



PA(r,*) / + -/_ 



eft 



PA(r,*) 



/+-/- 



(3.4) 



Here and below the expression PA means scalar production. 
Further we will use dimensionless velocity C in the form 

C = — = — - — A(r,t) = P- — A(r,£), 
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where P = — is the dimensionless momentum. 
Pt 

In linear approximation it is possible to replace the function / in Wigner 
- Vlasov integral by the absolute Fermi - - Dirac distribution, i.e. we put 
/ = If(P), where 

MP) = l + exp(P2-a) ' a = C ° nSt • 
Here Wigner — Vlasov integral (3.4) has the following form: 



= ^PAM) 



where 



If = fP(P) = 



If If 
1 



1 + exp 



2p q 



and pr = tuvt is the thermal electron momentum, or, 

1 



Jf 1 + e^"< 



Here 



or 



\ Px T ~y) + (^tj + { PzT —) 



n = 



2 



2 

Pt 



The linearization of the Wigner equilibrium function (3.2) we will carry 
out in terms of vector potential A(r, t): 

0/(0) 



/ (0) = / 



(0) 



a=o ~*~ <9A 



AM), 



A=0 



or, in explicit form: 



f^ = f F (P)+g(P)—PA(r,t), 

CPt 



(3.5) 



9^) = 



,P 2 -a 



(l + e p2 -"f 
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Considering decomposition (3.5), we will search for Wigner's function in 
the form: 

/ = f F (P) + g(P)—PA(r,t) + g(P)(PA(r,t))h(P). (3.6) 

CPt 

We receive the following equation 



PA(r,f) g(P)(v - iu + ikv)h(P) = 
2ie g(P)(u - v T kP) + ^M)(/+ - f F ) 



= PA 



Cp T 

From this equation we find 



ch 



PA(r, t)g(P)h(P) = PA(r,t)— - ^ - g(P)+ 



levr 



PA(r,f)- 



cpT v — ico + ifykP' 



(3.7) 

c/i v ' v — ico + ZfJ^kP 
With the help of (3.6) and (3.7) we construct the full distribution function 



= / (0) + 



/ = / (0) + g(P)PAh(P) = 
2ie co - v T kP A + iev T PA(/ + - fp 

CJ9t z/ — + ZUykP 



c/i z/ — icj + ivxkP ' 



or 



/ = /(°) + PA 



2ie lot — kiP . _ iel 

-,g{P) + 



f f fp 



cpx 1 — iior + ik{P ch 1 — -zcjt + zkiP 



(3.8) 



Here ki = k/, / is the electron mean free path, / = vtt, ki is the 
dimensionless wave vector. 

4. DENSITY OF ELECTRIC CURRENT 



We consider the connection between electric field and potentials 



c dt 



dr 
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or 



ILU 



E(r,*) = — A(r,*). 

c 

Hence, the current density is connected with vector potential as: 

j(r,t) = crt r — A(r,t). 
c 

By definition, the current density is equal to 



j(r,f) = e J v/- 



(2tt^) 3 ' 

Let's note, that the current density in the equilibrium state is equal to 
zero: 



j (0) (r^) = e /v(r,v,t)/(°)^ 



2p 3 d 3 P 



= 0. 



Indeed, considering that mean electron velocity in the equilibrium state 
is equal to zero, according to (2.4) we have: 

j (0) M) =eN^°\r,t) =0. 

Hence, with the use of equality (3.8) we have the following equality: 

j(r '' )= (2^/( PA ) V(r ' V '' )X 

d 3 P. 



x 



2ie ujt — kiP . _ iel 

g{P) + 



/f /f 



cpr 1 — ^wr + ikiP chl — iujt + ikiP 

Substituting obvious expression for the velocity into this equality 

p eA(r,f) p T P eA(r,f) 
v(r,r) = = , 

m mc m mc 

and, after linearization of it by vector field, we receive 

,4 



*£sfc/( A <'-*> p ) p '< 



X 



2ie cjt — kiP 
cpr 1 — ^wr + ik{P 



g(P) + 



id f+(P)-f F (P) 
chl — iujt + ik\¥ 



d d P. 



(4.1; 
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We notice that 

P =f J = P T — Pk + k 2 = 

^9 KV / \ 2 , 9 

= p T 2e; Pkl + fe) fc? - 

Let's copy the previous equality in the form 

j(r, *) = j^-rf- / [PA(r, *)]PS(P, Pk^P, (4.2) 

where 

cj9y 1 — zcjt + 2KiP ml - zcjt + zkiP 

A 

We take the unit vector ei = — , direct lengthwise the vector A. Then 

A 

the equality (4.2) we may write in the form: 

= 2i^M) I (Pei)PS( p PUl) rf3 Pi (4 3) 

In view of the symmetry the value of integral will not change, if the vector 
ei is replaced by any other unit vector e 2 , perpendicular to the vector k 1; 
i.e. 

A xki A xki 
62 = |A x ki| = Ah ' 
and A x ki is the vector product. 

Let's spread out a vector P in three orthogonal directions ei, e2 and n = 

ki 

ki 

P = (Pn)n + (Pei)ei + (Pe 2 )e 2 . 

By means of this decomposition we receive that 

(PA)P = A(Pei)P = 
= A(Pei)(Pn)n + A(Pei) 2 ei + ^(P ei )(Pe 2 )e 2 . 
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Substituting this decomposition in (4.3), and, considering, that integrals from 
odd functions on a symmetric interval are equal to zero, we receive that 

2ie 2 ^A(r, t) 

or, in the explicit form, 

2e/4A(r,f) 



y (Pe 1 ) 2 S'(P, Pn)d 3 P, (4.4) 



jM) = 



(27rh) 3 m 



eiP 



2ie wr — kiP 
cj9y 1 — icjr + ik.\P 



9(P)+ 



id f+(P) - fp(P) 



d 3 P. 



ch 1 — iwr + ikiP 

Here ei = is the unit vector directed lengthwise A. Therefore in 

view of symmetry 



2 r 2 

eiP) [S]d 3 P = / (e 2 p) [S]d 3 P = 



1 

2 



[S]d 3 P. 



eiPj + (e 2 P^ 

We will notice that the square of length of a vector P is equal 

P 2 = (P ei ) 2 + (Pe 2 ) 2 + (Pn) 2 , 

therefore 

(P ei ) 2 + (Pe 2 ) 2 = P 2 - = P 2 - (Pn) 2 = Pi 

where Pj_ is the projection of the vector P to direct, perpendicular planes 

(ei, e 2 ), and the vector n is the unit vector directed along the vector ki, n = 
ki k 

k\ k 

Hence for the current density we receive the following expression 



jM) = 



ep^A(r, t) 
(27rh) 3 m 



P 2 - (Pn)' 



2ie lot — kxP 



cpT 1 — icoT + ikiP 



g{P)+ 



iel f+(P)-fp(P) 
ch 1 — iujT + ikiP 



d 6 P. 
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Replacing the current density in the left-hand side of this equality by the 
expression in terms of field, we receive: 

ep^A(r, t) 



P 2 - (Pn)' 



x 



x 



2ie lot — kiP 
cpt 1 — iior + ik{P 



9(P) + 



iel f+(P)-f F (P 
ck 1 — iujT + ik\V 



d d P. 



5. ELECTRIC CONDUCTIVITY AND DIELECTRIC 

FUNCTION 

From the last formula we receive the following expression for the transversal 
electric conductivity in quantum non-degenerate plasma: 



2 4 



(27rh) 3 muj 



2 (Pki) 2 ^ 



P z - 



kl 



2 lot- kiP 



Pt 1 — ^wr + ikiP 



//+(P)-/-( P ) 



(5.1) 



ft 1 — icjr + ikiP 

We will transform expression for transversal conductivity and we will bring 
it to the form: 

2e 2 p\ 



(27rh) 3 mu J 
where &t is the thermal electron energy. 



PlSP 



1 — iut + ik{P 



With the use of the equality (2.3) we will present the previous formula in 
the form: 



47T/ 2 (a) 



[l-( W r)- 1 k 1 pl^(P)+ 
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/+(P)-/-(P) 



(P 2 - (Pn) 2 ) d 3 P 
1 — iujt + ik\V 



(5.1') 



Here the function /2(ck) has been entered above and in the absence of the 
scalar potential it is defined by equality: 



oo 



oo 



f 2 (a) = J x 2 f F (x)dx = J i * = \f ln(l + e°-*)dx. 



The quantity ctq is defined by classical expression for the static electric 
conductivity 



c^o 



Dielectric function we will find according to the formula: 

£tr = 1 H C<r- 

Substituting electric conductivity (3.1) into this equality we receive the 
expression for dielectric permittivity in quantum non-degenerate collision 
plasma: 

,2 • / 



e tr = 1 + HJ 



cj 2 47r/2(a) 



LOT - kiP 



^4- 



#(P)-/f(P) 



P 2 d 3 P 



(5.2) 



1 — iwr + ikiP 

We investigate some special cases of electroconductivity. In the long- wave 
limit (when k — > 0) from (3.1) we receive the well known classical expression: 

a tr (k = 0) = cr — = : — . 

V — 10J 1 — IUJT 

Let's consider the quantum mechanical limit of the conductivity in the 
case of arbitrary values of wave number, i.e. conductivity limit in the case, 
when Planck's constant h — > 0, and the quantity k is arbitrary. 

Now we consider the case, when values of the wave number are arbitrary, 
but Planck's constant converges to zero: h — > 0. 
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When the values of k are small we have: 

hk 



f±(P) = f F (P)±g(P)2P 



hence 



ff(P) - /f(P) = 2p(P)2P 



2mvT 
hk 



2mvT 



Therefore 



(uj - v T kP)g(P) + jeL[/+( P ) - /-(P)] = Uff( p). 

Thus, in linear approximation at small H (independently of the quantity 
k) for transverse conductivity we receive: 



_ -.classic 

°tr — &tr ' 



where 



classic = _^0_ f g(P)Pjd*P 
tT 47Tf 2 {a)J 1-lUT + lk!?' { J 

The expression (5.3) accurately coincides with the expression of the trans- 
versal conductivity for classical plasma with arbitrary temperature. 

Let's return to the expression (5.1'). We present it in the form of the sum 
of two components 

_ classic | quant fK A\ 

°tr = o tr + a tr , (5.4) 

where cr classic is defined by the equality (5.3), and second component af r uant 
corresponds to quantum properties of the plasma under consideration 



quant <^0 
< - 



47rf 2 (a) 



LOT 



£t [/^(p)-/f(p 



p2 j3d 

(5.5) 



1 — icoT + ik{P 



hw 

The quantum summand of r uant we will present in the form, proportional 
to a square of the Planck's constant h. 
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For this aim we use cubic expansion of af r uant by powers of h. We will 
remind, that in linear approximation by h, as it was already specified, the 
quantity <7^ uailt disappears. We will direct an axis x along the wave vector k. 

Let's expand the Fermi — Dirac distribution by degrees of dimensionless 
wave number q = - — = 1 2 , where kr = is the thermal wave number. 



We receive that 



hp mv c T 



h 



1 



/±(P) = f F (P) ± g{P)P x q - [g'AP)P z x + ^9(P) 



q 



± 



± 



P 



x „ + 



6 



Here g> p2 (P) = g'(P% g'UAP) = 9"(P 2 ) 

n P 2 -a 



g'(P 2 ) = g(P) 



1 



l + e P 2 -a 

Now it is easy to find the difference 



</'(P 2 ) = g(P) 



1 + e 



3 



P 2 -a 



UP) 



/+(P) - fp (P) = 2g(P)P x q + [g"(P 2 )P 2 x + -g\P 2 ) 
By means of this expression we find, that 



P x ^ + 



h P F 7^3 fe2 3 

^p) + |w(p)-ij(p)]=c(P)Sa 



+ 



where 



G(P) = P x [g"(P 2 )P 2 + -g'(P 2 ) 



3 



Substituting this expression into (5.5), we obtain, that the quantum summand 
is proportional to the square of Planck's constant and it is defined by expression 

kfu 3 f G(P)(P 2 - P%)d z P 



quant f^(j 

tr 247Tum 2 v^f2(a) J 1 — iur + ik\P x 

Let's similarly (5.4) we present the formula for calculation of dielectric 
permeability in the form 

classic 



1 i /classic , _ quant \ 

u~\ ) 



'tr 



i_ quant 

t tr 5 



(5.6) 
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where 



classic _ i i classic quant _ ^ quant 

t tr — i -|- o tr , t fr — a fr 

UJ UJ 



Thus, in an explicit form dielectric permeability of classical non-degenerate 
collisional plasmas is equal 

.classic = t , -<4 f g{P)PjSP 

tr + uj 2 ' 47r/ 2 (a) J 1-iwr + ikiP" { } 

And the component of dielectric permeability answering to quantum properties 
of non-degenerate collisional plasmas is defined by equality 

quant = 1 f -k 1 Pg(P) + (£, T /hv)(f+-f F ) 

£tr % uj 2 47Tf 2 {a)J 1-iuT + ikiP • [ } 

6. CALCULATION OF ELECTRIC CONDUCTIVITY AND 
DIELECTRIC PERMEABILITY 

In the expressions for classical and quantum components of the conductivity 
we can simplify several integrals. 

We break the triple integral to external one-dimensional integration by 
the variable P x from — oo to +oo and internal double integration by plane 
orthogonal to the axis P x in the expression (5.3). The internal integration we 
carry out in polar coordinates. Here we obtain that 

P 2 = P X + Pl, d 3 P = dP x dP ± , dP ± = Pj_ dP ± dx, 

where P± is the polar radius, and x is the polar angle. 
Thus we receive, that 

oo oo 2ir 

classic ^0 f jp f f g(P)PldP±d X 

a tr = A r , \ dP x 



47r/2(«) J J J 1 — + ik\P x 

oo 



where 
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Internal double integral we calculate in polar coordinates 

oo 2tt 

J J g(P) P'[dP ± d X = 7r\n(l + e a - p "). (6.1) 



Hence, the expression for the classical component is simplified to one- 
dimensional integral 

oo 2 
.lassie ^0 f H 1 + e a - P ^)dP x 

tr ~~4/ 2 (a) J 1 - icur + ihP x ■ 1 • } 



— oo 



The quantum item (5.5) we present in the form of the sum of two items 

quant . / r \ 

&tr =(Jl+(J2. (6.3) 

Here 

<jph r p x (p 2 - Pl)g{P) d 3 P 

47rf2(a)uT J 1 — iur + ik\P x 

and 

(T £ T f /f (P) ~ f F (P) (p2 p2 W 3p 

With the help of the equality (6.1) the expression for o\ can be rewritten 
in the following form 

oo 

_ aoh f P x ln(l + e a - p *)dP x 
Gl 4f 2 (a)cjT J l-iuT + ikiP x ' [ ' } 

— oo 

After change of variable 

the difference of integrals from (6.4) will be transformed to one integral and 
we receive 

»<T fc? f f F (P)(P 2 -P?)d 3 P 
SwMapT J (l-iwT + ik 1 P x ) 2 + (k 2 1 hv/4S. T y { ' 
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In the same way, as well as during the derivation of the formula (6.1), 
double internal integral in (6.6) we reduce to the one-dimensional integral 

oo oo oo 2ir 

J J f F (P)[P 2 -P*]dP ± = J J f F (P)PldP ± d X = 

— oo — oo 

oo 



= 2tt J P ± ln(l + e a - p "- p ") dP A 



o 

Now the expression (6.6) can be written in the following form (replacing 
a variable of integration P± = p) 



oo oo 

_ iaokf f f pln(l + e a -P- p *)dpdP x 
a2 ~ 4f 2 (a)cor J J (1 - iujr + ihP x f + {k\fiv/^ T f [ ' j 

-oo 

Thus, it is possible to present the expression for transverse conductivity 
in the form of the sum of one-dimensional (6.2), (6.5) and two-dimensional 
(6.7) integrals 

oo 

a f [1 - {h/ujT)P x } ln(l + e a ~ p *)dP x 



4/2 (a) J 1 — «wr + z/ciPa; 

—00 

OO OO 9 9 

zctq^i f f pln(l + e a -P- p *)dpdP x 



4f 2 (a)ur J J (l-iujT + ik 1 P x ) 2 + (k 2 hv/48. T ) 2 ' 

-00 

In the expression (6.6) for o 2 the thriple integral can be reduced to one- 
dimensional integral. For this purpose in (6.5) we pass to integration in 
spherical coordinates and present this expression in the form 



00 

io§k\ 

02 = 



00 

J f F (P)P 4 J(P)dP, 



4f 2 (a)uT 


where 

1 

j( P ) f (1 - ^ 2 ) dp 



1 (1-iuT + ikiPp) 2 + (klhis/4£ T ) 2 ' 
1 
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Let's designate temporarily 

i • h -u p w 

a = 1 — zcjt, o = z/ciP, a = — — — , 

and rewrite the integral J(P) in the form: 

l 

f (1 — fi 2 )dfi 

U 



7 (a + 6^) 2 + d 2 ' 

-l 



After change of variable a + b\i = t this integral will be rewritten in the 
form 

a+b 

b* J t 2 + d 2 

a—b 

This integral equals to 

2 d 2 + b 2 -a 2 1 (a + 6-id)(a-6 + id) 

^=-TT7H 7^ TTT^m 



6 2 6 3 2id (a + b + id){a-b-id) 

a (a + b — id)(a + b + id) 
b 3 (a — b — id) (a — b + zrf) ' 



or 



_2_ d 2 + b 2 -a 2 a 2 -{b-id) 2 ; a_ (a + b) 2 + d 2 
6 2+ 2idb* U a 2 -{d+id) 2 + b 3 U {a-d) 2 + d 2 ' 

Considering designations for a, b, d, we receive 

2 



J(P) ee J(P;wr, /ci) = 



(1 - wr) 2 + (hP) 2 - (huk 2 /48, T ) 2 (1 - iujt) 2 + (foP - hvkl/AZ T f 



k\P*{hv/2Z T ) (1 - wjt) 2 + (hP + hukl/48, T ) 2 

1-iuJT (l-iujT + ihP) 2 + {hvk\IM T ) 2 
+l (hP) 3 n (1 - zwr - z/ciP) 2 + {nvk 2 l^ T Y ^ ' 

Thus, the expression of quantum transverse conductivity is defined by 
one-dimensional integral 



4/ 2 (o;) / 1 - WJT + z/CiP x 



-00 
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ia k 2 



oo 



4f 2 (a)ur 



f F (P)P A J(P)dP, 



where the function J(P) is defined by expression (6.8). 

So, for electric conductivity and dielectric function we have received following 
expressions 



Oil 



Uf 2 (a) 



oo 



■— oo 



oo 



1 T\yln(l + e a T )dr y f uJ(u)du 

q xJ T — z/q xJ 1 + e u2 ~ a 

o 



(6.9) 



and 

£tr = 1 + 



x; 



4f 2 (a)x 2 



oo 



(x — qr) ln(l + e a T ") f u 4 J(u)du 

-dr — q 



t - z/q 



1 + e 



11 — (X 



-oo 



(6.10) 



In formulas (6.9) and (6.10) following designations are accepted 



k uj + iv lu 
q=-j-, z = x + iy=- , x = - . 

Kt KtVt KtVt 



V 



LUr 



y = 



krVT 



, x p 



= _2_ _ (z/q) 2 + (q/2) 2 - u 2 ^ (u - q/2) 2 - (z/q) 2 



u 3 q 



(u + z/q) 2 -(q/2) 2 



z ^(u- z/q) 2 -(q/2) 2 



[1 - r 2 )dr 



vfiq (u + z/q) 2 ~(q/2) 2 J (ut - z/q) 2 - (q/2) 2 ' 

-l 

Besides these formulas (6.9) and (6.10) we can give and other formulas for 
electric conductivity and dielectric function 



oo 



Oil 



y_yr^ ln(l + e a T ) ^ 



— oo 



4ih(a) J \q x J t - z/q 
dr 



oo 



oo 



4if 2 (a)x J (t- z/q) 2 - (q/2) 2 



J p\n(l + e a - T2 - p2 )dp 



-oo 



and 



£tr = 1 + 



Xp X 

4f 2 (a)x 2 q 



2 

q \ ln(l + e a ~ T )dr 
1 t I 



x J t — z/q 



-oo 
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oo oo 

xi f dr 



4/ 2 (^ J (r - z/ g y - (q ,2Y j " Ml + 

-oo 

7. THE SUM RULES 

Let's check up performance of one of the parities, named the rule /-sums 
(see, for example, [23]) for transversal dielectric permeability (6.14). This 
rule is expressed by the formula (4.200) of the monography 

00 



J Strict, v)(jjduJ = TTUJp. (7.L 



As shown in [22], for the proof of the parity (7.1) it is enough to prove 
performance of the limiting parity 

2 

£ tr iq,u, v) = 1 2 +0 (^)' w ^ OO. (7.2) 

From expression (6.10) it is visible that 

quant /' A 

By means of this limiting parity we will present the parity (6.10) at big u 
(u ^> 1) in the form 

tr + 47rf 2 (a)u 2 J 1-IUT + ik x P ' [ } 

From comparison (7.2) and (7.3) it is visible, that now it is required to 
prove equality 

iur f g{P)Pj d 3 P 
4irf 2 (a) J l-iujT + ikiP ~~ ~ ' 
or, considering once again, that u ^> 1, is required to prove equality 

-— L- f g{P)Pl d s P = 1. (7.4) 

Passing in equality (7.4) to spherical coordinates and integrating once in 
parts, we receive 

oo „ 

2 f e p - a P A dP 



3/ 2 (a)7 (l + e^ 2 -)2 



1. (7.5) 



o 
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Of justice of equality (7.5) we are convinced, having integrated once by 
parts. 

8. THE ANALYSIS OF RESULTS AND THE CONCLUSION 

The analysis of graphics on fig. 1-7 shows, that at great values q the 
quantum conductivity does not decrease to zero (as classical conductivity), 
and tends for collisional plasma to the finite limit. Really, from the formula 
(6.13) follows, that 

lim a tr = i— = i—. (8.1) 

q^oo X CO 

On Figs 1-7 all curves 1 answer to classical plasma, and curves 2 to the 
quantum plasma. 

On Figs. 1-3 are presented dependence \(Jtr /^ol quantum and classical 
conductivity from quantity q = k/kx at the various values of the resulted 
chemical potential a: a = (fig. 1), a = 6 (fig. 2) and a = — 5 (fig. 3). 

The analysis of graphics and numerical calculations show at negative 
values of chemical potential weak dependence Icr^/dol on quantity and 
strong dependence at positive values of chemical potential. From fig. 1-3 
it is visible, that while the module of classical conductivity has only one 
maximum, the module of quantum conductivity has also a minimum. The 
range of values the module quantum conductivity decreases at decrease chemical 
potential. 

On fig. 4 and 5 dependence Re(a tr / &o) is presented for quantum and 
classical conductivity. From the parity (8.1) it is visible, that the real part of 
quantum conductivity tends to zero at q — > oo the same as also the real part 
of classical conductivity. However, decrease of quantum conductivity occurs 
much faster, than decrease of classical conductivity that is visible from fig. 
4 and 5. At the real parts of quantum conductivity there is a maximum, the 
minimum is absent. 

On fig. 6 and 7 dependence Im (cr^/Vo) is presented quantum and classical 
conductivity on quantity q. Qualitatively imaginary part behaves the same 
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as also the module. From fig. 1 - 7 it is visible, that at q — > oo quantum and 
classical conductivity coincide. 

On fig. 8 dependence graphics |o"t r /cro| are presented quantum plasma 
from quantities q at various values of quantity x (dimensionless frequency 
oscillations of vector potential). The analysis of graphics shows, that with 
growth quantity x the range of values of the module of quantum plasma 
decreases, thus the maximum and module minimum are levelled, i.e. values 
of the module in points maximum and minimum approach, reducing range 
of values of the module of conductivity. 

In the present work the correct formula for calculation of transversal 
electric conductivity in the quantum non-degenerate collisinal plasma is deduced. 
For this purpose the Wigner— Vlasov— Boltzmann kinetic equation with collisional 
integral in the form of BGK-model (Bhatnagar, Gross and Krook) in coordinate 
space is used. 



Phc. 



1: Case: x = l,y = 0.01, a = 0. Dependence |o"t r /o"o| on quantity q. 




Phc. 2: Case: x = l,y = 0.01, a = 6. Dependence \crt r /cro\ on quantity q. 




Phc. 4: Case: x = l,y = 0.01, a = 0. Dependence Re(<%./cro) on quantity q. 




Phc. 6: Case: x = l,y = 0.01, a = 0. Dependence Im (atr/co) on quantity q. 



0.000 




Phc. 8: Case: x = l,y = 0.01, a = 0. Dependen |(o" tr /o"o)| on quantity q. 



34 



REFERENCES 

1. Klimontovich Y. and Silin V.P. The Spectra of Systems of Interacting 
Particles // JETF (Journal Experimental Theoreticheskoi Fiziki), 23, 
151 (1952). 

2. Lindhard J. On the properties of a gas of charged particles / / Kongelige 
Danske Videnskabernes Selskab, Matematisk-Fysiske Meddelelser. V. 28, 
m (1954), 1-57. 

3. Von Roos 0. Boltzmann — Vlasov Equation for a Quantum Plasma / / 
Phys. Rev. 119. JM (1960), 1174-1179. 

4. Kliewer K.L., Fuchs R. Lindhard Dielectric Functions with a Finite 
Electron Lifetime // Phys. Rev. 1969. V. 181. m. P. 552-558. 

5. Mermin N. D. Lindhard Dielectric Functions in the Relaxation-Time 
Approximation // Phys. Rev. B. 1970. V. 1, N>5. P. 2362-2363. 

6. Manfredi G. How to model quantum plasmas / / 
|ArXiv:quant-ph/0505004] 30 pp. 

7. Anderson D., Hall B., Lisak M., and Marklund M. Statistical effects in 
the multistream model for quantum plasmas// Phys. Rev. E 65 (2002), 
046417. 

8. De Andres P., Monreal R., and Flores F. Relaxation-time effects in 
the transverse dielectric function and the electromagnetic properties of 
metallic surfaces and small particles / Phys. Rev. B. 1986. Vol. 34,N 2 10, 
7365-7366. 

9. Shukla P.K. and Eliasson B. Nonlinear aspects of quantum plasma 
physics // Uspekhy Fiz. Nauk, 53(1) 2010; [V. 180. No. 1, 55-82 (2010) 
(in Russian)]. 



35 



10. Eliasson B. and Shukla P.K. Dispersion properties of electrostatic 
oscillations in quantum plasmas / / larXiv:0911 .4594V 1 [physics. plasm-ph] 
24 Nov 2009, 9 pp. 

11. Bhatnagar P.L., Gross E.P., and Krook M. A model for collision 
processes in gases. I. Small amplitude processes in charged and neutral 
one-component systems // Phys. Rev. 94 (1954), 511-525. 

12. Opher M., Morales G.J., Leboeuf J.N. Krook collisional models of the 
kinetic susceptibility of plasmas / / Phys. Rev. E. V. 66, 016407, 2002. 

13. Gelder van, A. P. Quantum Corrections in the Theory of the Anomalous 
Skin Effect // Phys. Rev. 1969. Vol. 187. m. P. 833-842. 

14. Fuchs R., Kliewer K.L. Surface plasmon in a semi-infinite free-electron 
gas // Phys. Rev. B. 1971. V. 3. N>7. P. 2270-2278. 

15. Fuchs R., Kliewer K.L. Optical properties of an electron gas: further 
studies of a nonlocal description / / Phys. Rev. 1969. V. 185. N 2 3. P. 
905-913. 

16. Dressel M., Griiner G. Electrodynamics of Solids. Optical Properties of 
Electrons in Matter / Cambridge. Univ. Press. 2003. 487 p. 

17. Wierling A. Interpolation between local field corrections and the 
Drude model by a generalized Mermin approach / / larXiv:08"T2 .3835vl 
[physics. plasm-ph] 19 Dec 2008. 

18. Brodin G., Marklund M., Manfredi G. Quantum Plasma Effects in the 
Classical Regime // Phys. Rev. Letters. 100, (2008). P. 175001-1 
175001-4. 

19. Manfredi G. and Haas F. Self-consistent fluid model for a quantum 
electron gas // Phys. Rev. B 64 (2001), 075316. 

20. Wigner E.P. On the quantum correction for thermodynamic equilibrium 
// Phys. Rev. 40 (1932), 749-759. 



36 



21. Tatarskii V.I. The Wigner representation of quantum mechanics // 
Uspekhy Fiz. Nauk. 26 (1983), 311-327; [Usp. Fis. Nauk. 139 (1983), 
587 (in Russian)]. 

22. Hillery M., O'Connell R.F., Scully M.O., and Wigner E. P. Distribution 
functions in physics: Fundamentals // Phys. Rev. 106 (1984), 121-167. 

23. Pains D, Nozier P. The theory of quantum liquids/ / Benjamin, New 
York, 1969. 



